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Figure 1: LO diagrams for bg — > tH 



1. Introduction 

Many models for new physics involve a more complicated Higgs sector than in the Standard 
Model. In the Minimal Supersymmetric Standard Model (MSSM) and other two-Higgs-doublet 
models (2HDM), one Higgs doublet gives mass to the up-type fermions and the other to the down- 
type fermions, with tanj3 = V2/V1 the ratio of the vacuum expectation values for the two doublets. 
The five physical Higgs particles in the MSSM include a light scalar, h°, a heavy scalar, H , a 
pseudoscalar, A , and two charged Higgs bosons, H + and H . A future discovery of a charged 
Higgs boson would constitute a definite sign of new physics. The Large Hadron Collider (LHC) is 
well positioned for a discovery of a charged Higgs [[]]]. 

A lot of work on higher-order QCD and SUSY corrections to charged Higgs production has 
been performed over the last several years, including next-to-leading order (NLO) calculations for 
bg -» tH |0, |, |], bb -> H+W~ [@, §, bb -» H+H §, qq -> H+H [§], as well as results 
for higher-order soft-gluon corrections for bg —* tH [^, |io| ]. 

2. Associated // and top quark production 

We begin with the dominant process at the LHC, which is associated charged Higgs and top 
quark production. The leading order (LO) process is bg — » tH and the corresponding Feynman 
diagrams are shown in Fig. [I[ The LO cross section is proportional to a a s (rattan 2 j8 + m 2 cot 2 /3) 
where nib is the bottom quark mass and m t is the top quark mass. 



Yukawa and SUSY electroweak corrections for this process were calculated in Qll| ] and 1 



loop SUSY corrections in Q12| , |13j|. The complete NLO QCD corrections were calculated in Ref. 
[§, [|, |J]. The QCD corrections were shown to be substantial, contributing up to 85% enhancement 
of the lowest order cross section and to reduce the scale dependence of the cross section. The 
NLO SUSY-QCD corrections are smaller in comparison, with their precise value depending on 
MSSM parameters [|, §]. 

To calculate the NLO QCD corrections, we have to include the one-loop virtual corrections to 
bg —>■ tH and also the processes with one additional parton: 

bg — > tH~g gg — > tH b qq — > tH~b bq — > tH q bq — > tH q 

bb -> tH-b bb tH-b 
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Issues with the calculation include the treatment of the bottom parton distribution, with a 
gluon splitting to bb in the collinear approximation, valid for small fr-quark pj. The diagrams for 
the process gg — > btH with a gluon splitting to bb are shown in Fig. ^. 

Work on matching the processes bg — » tH and gg — > btH [|14|] has been performed in [15, 



16, 17]. The use of matrix elements at large pj and parton showers at small pj, results in double 



counting for small pj when doing a simple addition. The matching performed in Q17Q involves an 
analytic double-counting subtraction term, and it can be implemented in event generators to give 
smooth differential distributions for all ft-quark pj. 



2.1 



tH near threshold 



Higher-order corrections to charged Higgs production near threshold have been calculated 
at next-to-next-to-leading order (NNLO) in Ref. [^] and at next-to-next-to-next-to-leading order 
(NNNLO) in Ref. [jT^]. 

For the process b(pi,) + g(p g ) — ► t{p t ) +H~{ph) we define the kinematical invariants s = 
(Pb + Pg) 2 , t = (pb— Pt) 2 , u = (p g — Pt) 2 , and S4 = s + t + u — mj — mji 2 , where m t is the top quark 
mass and m# is the charged Higgs mass. Note that, while we use the MS bottom quark mass in 
the coupling, we set m\, = in the kinematics. At threshold s<\ — > 0, and the soft-gluon corrections 
take the form [In 1 (s^/mfy/s^^ and can be resummed. For the order a" corrections / < In — 1. 
The leading logarithms (LL) are with I = 2n — 1 and the next-to-leading logarithms (NLL) are with 
l=2n-2. 

Near threshold these soft-gluon corrections are dominant and provide good approximations 
to the complete QCD corrections. The NLO and NNLO soft-gluon corrections were calculated at 
NLL accuracy in [§]. Furthermore, the NNNLO soft NLL corrections were presented in [10]. 

The calculation of these corrections is derived from the fixed-order expansion of the resummed 
cross section. Resummation follows from factorization properties of the cross section and is per- 
formed in moment space. We can write the resummed cross section as [^, [10[ |I8[ ] 



Cb%H-(X) = exp 



exp 



i JflF M 



exp 



£2 p (*00) 



r rVs/N fin 

H b g ^tH- (cf^^)) tfr-HR- faiy/s/fi)) exp / ^ 2Rer 



(2.1) 
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where the factorization scale is denoted by pLp and the renormalization scale by }JLr, and N is the 
moment variable. In the numerical results later we will set these two scales equal to each other and 



denote them by jit. The first exponent in Eq. ( |2.1| ) is 

dz- 



Y,E i {N i ) = -Y J C i f l dz 



z 



1 



1 



(l-z) 2 A 



dX a s (Xs) a,j((l 



-zfs) 



n 



K 



+ 0{at) (2.2) 



with Q = C F = (N 2 - 1)/{2N C ) for quarks and Q = C A = N c for gluons. 

The second exponent in Eq. ( |Q| ) involves the moment-space anomalous dimension y t n of the 
MS parton density, and the third exponent involves the QCD /3 function. H bg ^ tH and s bg ^ tH 
stand respectively for the hard-scattering function and the soft-gluon function, j^ 8 ^' 11 is the soft 
anomalous dimension, and its explicit form at one loop for this process is 



rbg^tH- 
1 s 



K 



C F ln 



-t + m, 
m t yfs 



-u + mf 
-t + mf 



in) 



+ ^(0 



(2.3) 



We then expand the moment-space expression of Eq. 1[ ) for the resummed cross section 
through NNNLO and invert back to momentum space. 
The NLO soft gluon corrections can be written as 



dt du 



K 



C3 



\n(s A /m 2 H ) 



5 4 



+ C 2 



(2.4) 



with F B the Born term, C3 = 2(Cp + Ca), and expressions for the other coefficients as given in 



The NNLO soft gluon corrections are 



daW{ 



54 1 



dt du 



-ci 



In 3 {s A /m 2 H ) 
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2 C3C2 "T C3 



ln 2 ^/^!)" 1 



S4 



+ ■ 



(2.5) 



where explicit expressions for subleading terms can be found in 
The NNNLO soft gluon corrections are [|l0|] 



dt du 



71° 





~\n 5 {s A /m 2 H )~ 
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5 2 



ln 4 (j4/m^) n 



+ ■ 



(2.6) 



and explicit expressions for the subleading terms and further details are given in JlOj]. 

2.2 H production via bg -> ?// at the LHC 

We now provide some numerical results for charged Higgs production in association with a 
top quark at the LHC. 

We first show that the soft-gluon corrections are dominant by comparing the NLO-NLL results 
with the exact NLO results that have been derived in Ref. [^J. We note that in [^] the reference 
scale chosen was pt = mn + m t . In most of the results below we will choose pt = ra#, which is a 
natural choice in our calculation. A cross section known to all orders does not depend on the scale. 
However a finite-order cross section does depend on the scale, though the dependence decreases 
as we move from LO to NLO, NNLO, NNNLO and so on. The work in [Q, ||] indeed showed 
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bg -> tH" at LHC S "=14 TeV n=m H -+m t 
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Figure 3: NLO exact and approximate K factors for H production at the LHC. 
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Figure 4: The total cross section (left) and K factors (right) for H production at the LHC. 



a reduction of scale dependence when the NLO corrections are added relative to the LO cross 
section. In fact, as we will see below, the higher-order threshold corrections further decrease the 
scale dependence, thus resulting in more stable predictions. But to make the comparison to [Q] 
we use a scale choice pL = m# + m t in Figure || and plot the K factors for H production at the 
LHC. The NLO-exact / LO curve shows the enhancement from the complete NLO corrections [§] 
while the NLO-NLL / LO curve shows the contribution of the NLL soft-gluon corrections at NLO. 
The two curves are close to each other and this is more easily seen from their ratio. The fact that 
the NLO-NLL / NLO-exact curve is very close to 1 (only a few percent difference) shows that the 
NLO-NLL cross section is a remarkably good approximation to the exact NLO result. 

In Figure ||, on the left, we plot the cross section versus charged Higgs mass for pp collisions 
at the LHC with \/S = 14 TeV using the MRST2002 approximate NNLO partem distributions 



functions [19] with a three-loop evaluation of a s . We show results for the LO, NLO-NLL, and 
NNLO-NLL, and NNNLO-NLL cross sections, all with a choice of scale jJ. = m# and with tan j8 = 
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K factors 


m H (GeV) 


NNLO-NLL 


NNNLO-NLL 


200 


1.34 


1.47 


300 


1.43 


1.53 


400 


1.49 


1.59 


500 


1.53 


1.65 


600 


1.57 


1.69 


700 


1.60 


1.72 


800 


1.63 


1.75 


900 


1.66 


1.79 


1000 


1.68 


1.81 



Table 1: The K factors for H production at the LHC. 



bg — > tH at LHC S =14TeV tanp=30 \tta g /2,2m g bg — > tH at LHC S =14TeV tan|3=30 m H =500 GeV 




m H (GeV) n / m H 



Figure 5: The scale dependence of the H cross section. 

30. The cross section spans three orders of magnitude in the mass range 200 GeV < m# < 1000 
GeV. The higher-order threshold corrections are positive and provide a significant enhancement to 
the lowest-order result. The cross sections for the related process bg — ► TH + are exactly the same. 

The right plot of Figure |] shows the relative size of the corrections as K factors at jj, = ra# . 
The NLO-NLL / LO curve shows that the NLO-NLL soft corrections enhance the LO result by 
25% to 48% depending on the charged Higgs mass. The K factors increase with higher masses, 
as expected, since then we get closer to threshold. With the NNLO-NLL corrections added we get 
an enhancement over the LO result of 34% to 68%. Adding further the NNNLO-NLL corrections 
provides an enhancement ranging from 47% to 81% over the LO result. 

In Table 1 we show the explicit numbers for the NNLO-NLL and NNNLO-NLL K factors for 
specific values of the charged Higgs mass with pL = wir- 

We next study the scale dependence of the cross section. On the left plot of Figure |5| we show 
the cross section at the LHC with tan /3 = 30 as a function of charged Higgs mass with two different 
choices of scale, /I = mn/2 and 2m#. The scale variation of the LO cross section is quite large. 
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The variation at NLO-NLL is smaller, and at NNLO-NLL and NNNLO-NLL it is very small. In 
fact the two NNLO-NLL curves are on top of each other for much of the range in as are the 
two NNNLO-NLL curves. 

On the right plot of Figure g, we show the scale dependence of the cross section for m# = 500 
GeV and tan/3 = 30 over a large range in scale, 0.1 < \ljmn < 10. The higher-order threshold 
corrections greatly decrease the scale dependence of the cross section. The NNNLO-NLL curve is 
relatively flat. This can also be demonstrated by calculating at each order the ratio of the maximum 
value to the minimal value of the cross section over the /J. range. We find 

a max /a rain = 3.39 1.50 1.38 1.32 

T T T T 

LO NLO-NLL NNLO-NLL NNNLO-NLL 

We see that with progressing order the ratio decreases and gets closer to one. 

In Figure [] we plot the dependence of the H cross section on tan/3 (left) and the top quark 
mass (right) with pL = rriH = 500 GeV. The tan /3 variation in the left plot is over the range 1 < 
tan/3 < 50, and the cross section is at a minimum near tan/3 = 8. We note that the tan/3 depen- 
dence arises in the factor m\ tan 2 j8 + mf cot 2 /3 , and the tan j8 shape is the same for all curves. The 
dependence on tan /3 is large, spanning two orders of magnitude in the range shown. The depen- 
dence of the cross section on the top quark mass is shown in the right plot for tan /3 = 30. For 
heavier top quark masses the cross section decreases. We see that the dependence is mild so that 
the current small experimental uncertainties on the top quark mass do not play a major role in the 
total uncertainty of the charged Higgs production cross section. 



3. Other charged Higgs production channels 



We briefly discuss some other production channels for charged Higgs production, including 
associated production with a W boson, and pair production. 
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Figure 7: LO diagrams for bb — > H + W 
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Figure 8: LO diagrams for gg — > . 



3.1 Associated H + and W production 



Charged Higgs bosons can be produced in association with W bosons (see, e.g. [Q, ^, 20, 21 ]). 
The LO processes are gg — ► and M> — > . LO diagrams for bb — > are shown 

in Fig. [7| 

Complete NLO calculations for bb — > were presented in §]. In Ref. [21 ] the above 

production process followed by a leptonic decay of H + and a hadronic decay of W~ was studied 
as a possibility for observing charged Higgs bosons at the LHC. 

3.2 Charged Higgs pair production 

The LO processes for the production of a H + , H pair are gg — > H + H , bb — ► H + H , and 
— ► H + H with g a light quark. 

LO diagrams for — * H + H~ , involving loops with top and bottom quarks, are shown in 
Figure H. Studies have been made for charged Higgs pair production via this process at the LHC 
(see e.g. @ ^3[|§). 

LO diagrams for bb — > H + H are shown in Figure |9| The NLO corrections for this process 
have been presented in §]. The NLO corrections for the Drell-Yan type process involving light 
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Figure 9: LO diagrams for bb — > H + H 



quarks, qq — > H + H , were also presented in [§J. At LO the Drell-Yan type process proceeds via 
an s-channel 7 or Z boson, similar to the left diagram of Figure ^. 

The relative contribution of the gg, light qq, and bb channels to charged Higgs pair production 
at the LHC depends on the values of tan/3 and Higgs mass. At tan/3 = 10 the light qq conuibution 
is by far dominant but at tan jS = 50 the gg contribution dominates [Eh. 



Ref . [ 25 ] provides a study of the associated production of a charged Higgs pair with a bb pair, 
gg — > bbH + H , which is the dominant pair production mode at large tan/3 and relevant for the 
determination of triple-Higgs couplings. 

Further calculations of NLO and higher-order corrections will be crucial in reducing uncer- 
tainties in the theoretical predictions for cross sections and differential distributions for charged 
Higgs production. 
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